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Abstract 

We discuss a new relation between the low lying Schroedinger wave func- 
tion of a particle in a one-dimentional potential V and the solution of the 
corresponding Haniilton-Jacobi equation with —V as its potential. The 
function y is > 0, and can have several minina {V = 0). We assume the 
problem to be characterized by a small anhamornicity parameter g^^ and a 
much smaller quantum tunneling parameter e between these different min- 
ima. Expanding either the wave function or its energy as a formal double 
power series in g^^ and e, we show how the coefficients of g^"^e^ in such an 
expansion can be expressed in terms of definite integrals, with leading order 
term determined by the classical solution of the Hamilton-Jacobi equation. 
A detailed analysis is given for the particular example of quartic potential 
V = \g\x^ - a'f. 

PACS: ll.lO.Ef, 03.65.Ge 



1. INTRODUCTION 

We discuss a hitherto unexplored link between the low-lying Schroedinger wave function of 
a particle in a potential +V and the solution of the corresponding Hamilton-Jacobi equation 
with —V as its potential. In this paper, we restrict our study to the one dimensional problem. 
Extension to multi-dimensional space will be given in a separate publication. 

Let $(a;) be the solution of 

H^x) = E^x), (1.1) 

where 

//=-i^ + VW (1.2) 

describes the Hamiltonian of a non-relativistic particle of unit mass. Assume 

V{x) > 0; (1.3) 

however, its minimum V{x) = may occur at more than one point. 
Introduce a scale factor g'^ by writing 

V{x)=g%{x), (1.4) 

and consider the case of large g. We express 

$(a;) = e-^(") (1.5) 

in terms of a formal expansion: 

S(x) = gSoix) + Si(x) + g-'S2ix) + g'^S^ix) + ■■■ (1.6) 

and likewise its eigenvalue 

E = gE^ + E^+g-^E2 + g-^E^ + ---. (1.7) 



Substituting ( |1.5|) - (|1.7| ) into the Schroedinger equation (|l.lj ) and equating the coefficients 



of g " on both sides, we obtain 

S[? = 2v, 
SoSi = -So-Eo, (1.8) 

etc., where S^ = dSn/dx and S" = (fSn/dx^. For each n > 1, the (n + l)th line of (|0|) 
determines S'„ uniquely in terms of the S^ for m < n. But the first line determines Sq only 
up to an ambiguity in sign. By continuity, the ambiguity makes itself felt only at a global 
maximum of —v{x), that is where v(x) vanishes quadratically and we must decide whether Sg 
should behave analytically, changing sign as it passes through the vanishing point, or whether 
S'o should have a kink and take the same sign on both sides of the vanishing point. 

If v(x) vanishes only once, say at x = a, our choice is clear. Sq should be positive for 
X > a and negative for x < a, so that So itself will increase on both sides as we go away from 
X = a, and the wavefunction in ( |1.5| ) will vanish at infinity in both directions. Thus we should 
take the "analytic" square root in the first line of (|1.8| ). 

If there are M > 1 values of x where v(x) vanishes, let us call one of them {x = a), 
where we want the wavefunction to have its major peak, the primary maximum of —v, and 
the others remote maxima. At the primary maximum Sq should behave analytically as before. 
But if we allow Sq to change sign also at a remote maximum, we shall find So decreasing 
as we move farther on the other side, so that the wavefunction may grow to another peak 
competing with the primary one, or may blow up at infinity. 

Thus, to have a well-behaved wavefunction with no remote peaks, it appears that Sq should 
change sign only at the primary maximum, and should have a kink at the remote maxima. 

Further support for this choice arises from the generalization of our problem to many 
dimensions. The first line of ( |1.8|) will then read (ySo)^ = 2f , and at the primary maximum 



it will be natural to make ySo a vanishing vector with positive divergence. Elsewhere Sq 
will be found by integrating this (Hamilton-Jacobi) equation along a family of trajectories 
radiating out from the primary maximum. In general (excluding a set of measure zero) these 
trajectories will not encounter any remote maximum. If we follow a trajectory that comes 
very near to a remote maximum, we shall see —v rising almost to zero but not quite. Thus, 
Sq (differentiation along the trajectory) will not change sign, but will not quite have a kink. 
Passing to the limit, we see that in order to make Sq continuous in all dimensions we should 
allow S'o to develop a kink on the singular trajectory that passes through the remote maximum. 

In the present paper we shall be concerned more with the quantum mechanical wavefunction 
than with the Hamilton-Jacobi equation, as we wish to develop corrections to all orders in g^^ 
and in e, the tunneling exponential. We shall find it convenient to introduce two types of 
wavefunction denoted respectively by </> and tp. Generically, $ refers to either (p or -ip. 

The groundstate energy will have an M-fold quasi-degeneracy with a spread of 0(e). 
Correspondingly, there are M eigenstates (called quasi-groundstates) of H. A 0-function 
will be an exact solution of the Schroedinger equation with an energy E that is not a true 
eigenvalue but lies within the groundstate fine structure. A 0-function cannot behave well 
at both +00 and —00. On the other hand, a ^-function is a linear combination of the true 
quasi-groundstate eigenfunctions. A ^-function behaves well at both infinities, but is not an 
exact solution of the Schroedinger equation for any energy (except when ip happens to be just 
one of the M quasi-groundstate wave functions). 

In this and the following section we shall study only the leading order in e, and for this 
purpose we can restrict ourselves to the region excluding any remote maximum, i.e., to x such 
that there is no second vanishing of v between x and a. In this region 0- and ^-functions 
are almost indistinguishable, and all the assertions in these sections will apply indifferently to 
both, except for the exact equations (P^), (|2.1C1| ) and ( |2.14|) involving the Wronskian. 



In section 3 ff., we shall have to go beyond the remote maxima, and for Af = 2 at least we 



shall find that for a 0-function the energy E can be tuned so that the Hamilton-Jacobi solution 
with analytic Sq remains a good approximation, whereas for a ^-function the coefficients of 
the eigenfunctions can be tuned so that the Hamilton-Jacobi solution with a kink in Sq at the 
remote maximum is good. 

It is convenient to write the first line of (|1.8| ) as 



i , dog ,2 



v{x) = 0+, (1.9) 



2 dx 

which denotes the corresponding Hamilton-Jacobi equation with 0+ as the total energy and 
—V the potential. For its solution, we follow Hamilton's action-integral by starting at the 
primary maximum of —v{x), say x = a; i.e., 

-v{a) = 0. (1.10) 

At this starting point, the magnitude of Sq is infinitesimal because the total energy is 0+. 
However, we will allow Sq to change sign at x = a in accordance with the discussion following 

( p..8|) . Thus, Sq is positive for x > a, and negative for x < a. In the limit OH > 0, Sq = 

at X = a and also at any other (degenerate) maximum of —v{x), say x = b and b ^ a. The 
limiting Sq becomes analytic at x = a, but has a kink at x = b. The resulting $ = e~^ 
describes then a ip-fur\ct\or\. 

Impose the "normalization" condition 

<l>(a) = 1 (1.11) 

and, correspondingly, 

So(a) = Si(a) = S2(a) = ■•• = 0. (1.12) 

The solution of the Hamilton-Jacobi equation with the sign choices described above is 

X 

So{x) = ±J[2v{y) + {0+)]-^dy (1.13) 



with [■ ■ ■]2 denoting the principal value, the "+" sign for a; > a and the "— " sign for x < a, 
so that 



Sa(x) > 



everywhere. 



As we shall see, the requirement that Si(a;), 82(0;), ■■■be regular at x = a determines 
Eq, El, E2, ■ ■ ■■ 

As a first illustration, consider the simple harmonic oscillator with 



V{x) = ^x^ 



There is only one minimum of V{x) at 



X = a = 0. 



1.14) 



Thus, (|r|) and ( p^ - (gJ^) lead to 



Si(x) = 82(0;) = ■■ ■ = 0, 
Ei=E2 = --- = 0, 



S{x) = gSo = -gx^ 



and 



;i.i5) 



E = gEo = -g. 



The solution of the Hamilton-Jacobi equation 

1 C?S 2_^ 2_ 

gives the exact groundstate Schroedinger wave function 



;i.i6) 



;i.i7) 



^(x) = e~^(^'^ = e"5S^'. 



:i.i8i 



While the above approach resembles the familiar WKB method, it is different. The latter 
is good for excited levels, but not for the groundstate: the WKB wave function has an 
unnecessary "classical turning point" at x = ±g^^ (even though it does give the correct 
eigenvalue). Note that WKB calls for introducing the energy E on the right hand side of the 
first line of (|1.8| ), whereas we defer it to the second. 

As a second example, we consider the quartic potential 

V{x) = ^g\x'-ay. (1.19) 

The minimum V{x) = is realized at two separate points 

X = ±a. (1.20) 

In this problem, in addition to the small anharmonicity parameter 

{ga^' « 1, (1.21) 

there is a much smaller associated parameter 

e = e-l'^"' (1.22) 

that characterizes the tunneling between the two minima x = a and —a. We shall start with 
the two solutions So(±) of the Hamilton-Jacobi equation ( p..9|) , 

^[^So(±)F-^(x^-a^)^ = 0+ (1.23) 

by requiring 

So(+) =0 at X = a 

and (1.24) 

So(— ) = at X = —a. 



Hence, 



X 



So(+) = ±/[(y2-a2)2 + (0+)]trfy 

a 

and (1.25) 



in accordance with (|1.13| ) and the sign choices described below ( p..lO|) , so that So(+) and 



So(— ) are positive everywhere. In addition, they are related to each other by 

X -^ -X, So(±) ^ So(t). (1.26) 

Let H be the Hamiltonian associated with the quartic potential ( p..l9| ): 

*-i^4«V-V. (1.27) 

and ^even, tpodd its even and odd ground-states with eigenvalues Eeven and Eodd- Write 

^ 'Zy^odd ~r -L-'even): 

A = -{Eodd — Eeven) (1.28) 



and 



2 



i^±{x) = -[lpeven{x) ± tpoddix)]. (1.29) 



Thus, 



ip+{x) = i;_{-x) (1.30) 

{H - S)ij^{x) = -Aij^ix) (1.31) 



and 



(^;^_ _ ^'_^^y = 2A ■ (^2 _ ^^)^ (1.32) 



where, as well as throughout the paper, the prime denotes ^. 

We may expand S and A both formally as the following double series: 

S = gaY: C^„(<7a3)-'"e-t'^3»' (1.33) 

m,n 

and 

A = A(-g^a'')h-l<^^' Y Cmniga'^y'^e-^^'''''' (1.34) 

TT — ' 

where Cmn and Cmn are numerical coefficients. The main body of this paper is to show how 
by using the solution ( |1.25|) of the Hamilton-Jacobi Equation as the zeroth approximation to 



the Schroedinger wave function ( pTsp , we can obtain explicit expressions for these coefficients. 
In Section 2, we derive 

Coo = Coo = 1, 

Cio = -\, C20 = -9/2^ (1.35) 

^30 = -89/2^ etc. 

The general formulas for Cmn and Cmn are obtained in Section 3 and 4. Discussions for other 
examples of V{x) are given in Section 5. In the Appendix, we analyze further the difference 
between a ^/^-function and a ^-function. 
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2. ZEROTH AND FIRST ORDERS IN e 

We continue our analysis of the quartic potential case: 

1 



V{x) = -g'{x'-ar. 



To the zeroth order in 





e = e 




-^even 


= ^odd 


= E. 


T3^), 





(2.1) 

We write in the notation of ( |1.33| ), 

00 
E = gaY. C^.iga')-^. (2.2) 

m=0 

(This definition is a bit fuzzy since the series is asymptotic. In the next section we shall define 
E precisely. For the purpose of this section, it suffices that E = ^{Eeven + E^dd) correct to 
0(e) and that ( |2.2|) follows as a consequence. We define 0-j- as the exact solution of ( |1.1| ) 
with 0+(cx3) = 0_(— 00) = 0. But since E is not an eigenvalue of H, 0+(— 00) and 0_(oo) 
are both divergent.) 

For clarity, we reserve ip± ^or the exact wave functions given by (|1.29| ), and use 

ip+{x) = (p+{x) for X > —a, 

^_(x)^0_(x) for x<a (2.3) 

in the approximation of neglecting e. 

Following ( |1.5|) and setting the generic $ as (f)± we write 

(f)^{x) = e-^(±) (2.4) 

with 
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0+(x) =0_(-x), 

0+(oo) = 0_(-oo) = (2.5) 

and 

0+(a) = 0_(— a) = 1. 



In accordance with (|1.6| ) - (|1.7| ) and writing S as 5'(±), we expand 



Si±)=gj:g-"'SU±) 

m 

and (2.6) 



■'m 
m 



where 

Sm{,+) = at a; = a 

Sm{-) = at a; = -a, (2.7) 

and the corresponding S'^{+) and 5'^(— ) are finite at a; = a and x = —a, respectively. As 
mentioned before (after ([2.2|)), we have 

{H - E)<P± = 0, (2.8) 

which leads to 

(0V0- - 0'_0+)' = (2.9) 

and therefore 

0'_^0_ — 0'_0+ = A = const. (2-10) 

By using the Hamiltonian ( p..27| ) and the set of equations ( p..8|) but replacing S^ by S'm(±), 
we derive 

12 



\x + a\ >0{-^) for Si- 



\x-a\ >0(^) for Si- 



Sq{+) = i(x-a)2(a; + 2a) 

S2i + ] 



1 

x+a 



x+2a 



16a^ 4a2(x+a)2 

a:+3a 



S2{ + ) '- 4a2(x+a)3 

C /'_|_\ _ 71 _ 9x^+36ax^+57a^x+40a3 
•JSl-^J — 29a6 64a5(x+a)4 

C/('_i_\ _ 9x^-\-45ax'^+99a'^x+103a^ 
'^3\^J — 64a5(x+a)5 



^4( + ) 



387 



Soi-) = lix + afi-x + 2a) 

S'oi-) = -ix' - a') 
SiH = In^ 



S2H = 
S'2{-) 



1 

x—a 



3 
16a^ 



-x+2a 
4a2(-x+a)2 

—x+3a 



Aa^{—x+a)^ 



'S'3( — ) — T^g-fi + 



29 a6 



9a:^ - 36ax2 +57a^x-40a^ 

64a5(-x+a)4 



S'{ 



9x'^-45ax^+99a^a:-103a^ 
64a5(-a'+a)5 



5^ + 1 



2iia9 

267x^+1602ax'^+i094:a^x^+57ma^x^+4779a^x+2068a^ 
3-2^a»{x+a)^ 

89x^+623ax^+1958a^x^+3594a^x^+4121a*x-\-25^3a^ 



512a»(x+ay 



Si{-) = Si{+){x -^ -x) 



S',{-) = -Si{+){x ^ -x) 



(2.11) 



and 



E = ga 



1 



9 



4a2 2^(2a)5 2^2(2a) 



89 ^/ 1 N 



A = 8^aV3»'^ [1 



53 



+ 



(2.12) 



:2.13) 



8ga^ 256(^a3)2 
This expression of E gives the coefficients Cio — C30 quoted in ( |1.35D . 

We note that for a; < —2a, ( |2.11|) gives a negative S'o(+), and as expected, when x -^ — cxd, 
'S'o(+) -^ —00 and therefore 0+(— 00) is divergent; likewise, as x ^ +00, So{—) -^ —00 and 
0_(oo) is divergent. 

In order to evaluate A, we return to the exact wave functions z/'i. given by ( p..29| ), and 
integrate ( |1.32| ) from —00 to x. Since ipevenix) and ipoddix) are eigenfunctions of H, both 
are well behaved at x = ±00, so are i'±{x) of ( p.. 291) . Hence, ( |1.32|) gives 
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ij'4x)Mx) - ^Ux)^+(x) = 2A / l^Uy) - ^l{y)]dy. 



(2.14) 



According to ( |2.3|) - ( |2.4|) and (|2.6| ), the zeroth order of the exponents '\n ijj± = e ^^"^^ are 



determined by ( |2.11J ) provided x > —a + 0{-^) for i>+{x) and x < a — 0(^=) for ^p 



[X 



Take 



a < X < 



(2.15) 



and also 



/ 1 ^i 

X + a » ( — )2. 
ga 



Hence, the left hand side of ( p. 141) is given by ( p.l3|) , while its right hand side is 2A times 



my)-i^Uy)]dy 



2S(^)dy + 0ie) 



which leads to 



and i.e. 



-2ga{y+a)2 



dy[l+Oi—)+0{e)] 
ga-^ 

;^)^[i+o(^)+o(6)], 

2ga ga-^ 



, ,2qa. 1 

TT 



A = 4^a2(^)^e-t^'^'[l + 0{\) + 0{e)]. 



(2.16) 



(2.17) 



Thus, Coo = 1 in ( |1.34| ). By using the first line of ( p.l6|) and expanding exp[—S 



g^^S2{—) — g^'^Ss{—) ■ • ■] as a formal series in g^^, we can evaluate cio, C20, ■ ■ ■ in terms of 
moments of Gaussian integrals. 

In the next two sections, we develop a systematic method which enables us to evaluate all 
other Crnn and Cm„ in (|1.33| ) - ( |1.34|) in terms of definite integrals. 
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3. DEFINITIONS AND PRELIMINARIES 

In this section we introduce several definitions and establish some useful formulas for the 
quartic potential case. 

In (|2.3|), the i'±{x) on its left hand side are defined in terms of the eigenfunction ipeven 
and 4'odd of the quartic-potential Hamiltonian H, in accordance with ( p..27| ) - ( |1.29|) , whereas 
the 0±(x) on its right hand side are introduced as solutions of 

{H-E)Mx) = 0, (3.1) 

with E between the two lowest eigenvalues -Eetjen and Eodd, 

Eodd > E > Eeven, (3.2) 

and (p± satisfying the additional conditions 

0+(x) =0_(-x) 

and (3.3) 

0+(oo) = </)_(— oo) = 0. 

It follows from the Sturm-Liouville Theorem that (p+{x) and </>-(a;) must have one (and only 
one) zero each at a finite x. Let 

0^(-«) =0_(a) = 0. (3.4) 

When E decreases from Eodd to Eeven, « moves from z = to oo. It is convenient to choose 
a to be near a, with 

\a-a\<0(^). (3.5) 

x/ga 



As we shall see in Section 4, (|3.5|) is equivalent to constraint 
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|^-^(^e.en + ^odd)|<0(e2) (3.6) 



where e is given by ( ^^1]) . Because of {q^, (0'+*^- — 0'_0+)' = is exact, valid to all orders 
ofe = e-y^\ and so is ( ^1^) : 

(t>'_^_(f)- — ((>'_((>+ = A = const. (3.7) 

Likewise, the asymptotic expansion ( |2.13| ) 



holds at all x. From (|]3|) - {^ and (|]^, we have 



/ -^jT for X < a 
0+(x) = A0_(x)-<{ ~»J-'" (3.8) 

and 

[ J^ for x>-a 

For |x + a| > 0{^=), the exponent S'(+) of 

0+(a;) = e"'^W (310) 

has an asymptotic expansion given by ( |2.6|) and ( [2.11] ). Through 

0^(-x) = 0_(x) = e-^(-), (3.11) 

we have the corresponding asymptotic expansion of S{—) for |a; — a| > 0{^=). These 
expressions and the Schroedinger equation ( p.l| ) lead to the previously derived asymptotic 
expansion (|2.12| ) for E, independent of a. To the zero*^ order in e, we note that for ^+(a;) of 
( p..29| ), the asymptotic series is valid only for x > —a + 0{-^); however, as will be discussed 

16 



later, for (p+{x) the same asymptotic series also holds for x < —a — 0{^=) if a satisfies 
( [3.5|) . A similar difference exists between ip^{x) and (p-{x) through x — > —x. 

Different choices of ±a give different 0±(x), and therefore also different E. The asymp- 
totic expansion for E implies that these different E become the same to the zeroth order in 
e. (It will be shown in Section 4 that actually the constraint ( [3.5|) determines E to the first 
order in e.) 

Similarly, for x > —a + 0{^=), the a-independent asymptotic expansion ( [2.11|) for S{+) 



implies that these different (p+{x) are also the same if we neglect the 0(e) correction. For 

x<-a + 0(^=), (3.12) 



{x) is most conveniently given by (^]8|). We note that for x < a — 0(^=), which includes 



go.' 



the above region ( |3.12|) , the asymptotic expansion of (|2.11J ) for S{—) determines 0-(x) to 



the zeroth order in e. Thus, in the region ( [3.12|) , under a ^ a + 5a, we have from (^]8| 



5^+{x) = \[(p4x)/(pl{-a)]Sa, (3.13) 

in which \S(f)- = O(e^) is neglected. Since [1/0^(|/)] is minimum (for y < a) at y = —a, the 
restriction ( |3.5| ) limits the corresponding 50+ to 0(e) or smaller. Likewise, through x -^ —x, 
the corresponding 5(f)- is also 0(e) or smaller. Furthermore, for a satisfying the restriction 
(p. 51), in the region x is less than —a by 0(^=) 



X < —a — 0( )■ 



through the upper expression of ( p.8|) one can use the asymptotic expansion of </)„ = e^"^*^"^ 
to derive the asymptotic expansion of 0+ = 6""^^+^ in accordance with ( p.ll|) . A similar 
statement holds in which x -^ —x and the roles of 0_ and 0+ are reversed. 

In Section 4 we shall relate the eigenfunctions ipeven and ijjodd to these 0-|-. For this purpose 
it is useful to introduce two Green's functions ("right" and "left"): 
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<x\GR\y>=j- I (3.14) 

for X > y 



2 for y > X 

<x\GL\y>=j-l (3.15) 

(j)+{x)(j)^{y) - (j)^{x)(j)+{y), for x > y. 



and two corresponding step functions: 



and 



Hence, 



— 1 if y > X 
SRix-y) = { (3.16) 

if X > y 



if y > X 
SLix-y) = { (3.17) 

1 if X > y 



2 

< x\GR(^L)\y >= j[(p+{x)(f)^{y) - (f)_{x)(f)+{y)] ■ Sr(l){x - y). (3.18) 



Both Sr and Sl satisfy 

d 

-^SR(^L){x-y)=6{x-y). (3.19) 

Since 

[0+(x)0_(t/) - (j)_{x)(p+{y)]6{x - y) = 
we find 

— <x\GR(^L)\y> = -[0'+(^)0-(l/) -0-(^)0+(l/)]'^R(L)(^-l/)> (3-20) 

1 (9^ 2 1 

~2 ar2 < ^I'^^wb > = -y{-i^WA'^)'^~^y) - <^'-{^)<^Av)\Sr{L){x - y) 

^\{-\)\<^\{x)<^.{y)-^'_{x)^^{ymx-y) (3.21) 

18 



and 

{H-E)< x\GR^L)\y >= -jW4x)<p4y) - <p'_{x)My)] ■ Si^ -y) = -K^ - y)- (3.22) 

Note that at a fixed x, as y — i> oo 



< AGrW > — ^ -j<f)+{x)(j)^{y) — > oo 



(3.23) 



and as y — > — oo, 



< xlGily > — > --(f)-{x)(f)+{y) — > oo 



(3.24) 



Because of (jT^) and { ^J^ - (|m| ), at all x we have 



^+[x)(p--[x} 



A 



si-y-si+y^ 



(3.25) 



furthermore, for 



X - a \ » 0( 



(3.26) 



we find, by using ( p.llj ). 



si-y - sM - Ma^ -V j^ - %if^ + 0(1) 



(3.27) 



Therefore, integrals such as 



< x\GR\y > (p+iy)dy 



(3.28) 



and 



< x\GL\y > (t>-{y)dy 



(3.29) 



are convergent. 

Next, we examine the equation 
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if0+ = -^0'| + r0+ = E0+, (3.30) 

in which E does not have to be an eigenvalue of H. Write, as before, 0+ = e"'^^"^-' but 
introduce a new function 9, defined to be 

e = S{+)'-V2V. (3.31) 

Thus, 9 satisfies 

-(9' -9^)-V2V9 + --^ = E (3.32) 

which, for the quartic potential 

becomes 

^9' -9^)-g{x^-ay + gx = E. (3.33) 

The boundary condition 0+(oo) = gives, for large x, 

X gx^ x^ gx^ 

The asymptotic behavior of 0+(a;) as x — >^ oo is therefore insensitive to E; the same holds 
for <\)-{x) as s ^ — oo. Consequently, the convergences of (p.28| ) - ( p.29| ) also imply the 
convergences of 



<x\GRi^V)\y> i)even{y)dy (3.35) 

— oo 

and 

oo 

<AGR{L)\v>i^odd{v)(i-V (3.36) 
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According to ( 3.14| ) - (|3.15| ), each of these Green's functions Gr and Gl consists of two 
terms, — 0+(a:;)0_(y) and 0_(x)0+(y). It is useful to examine their relative magnitudes. Since 
0+(— a) = 0-(a;) = 0, by using ( p.ll| ) we see that (j)+{x) > for x > —a and (t)-{x) > 
for X < a. Because 

r(j)-ix) , A 



0+(x) 0^(X) 

for any x < y we have 



< 0, (3.37) 



Hence for —a < x < y 

(t>-{x)(t>+{y)>(t>-{y)M^)- (3.39) 

For 

— a<x<y<a, (3.40) 

both sides of ( [3.39| ) are positive, which gives 

< 0-(x)0+(y) - 0-(i/)0+(x) < 0_(x)0+(i/). (3.41) 

For 

a <x <y, (3.42) 

both sides of ( [3.39| ) are negative. We have 

< -0_(x)0+(y) < -0+(x)0„(y) (3.43) 

and therefore 

< 0_(x)0+(i/) - 0+(x)0_(y) < -0+(x)0„(y). (3.44) 
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These inequalities help us to study the bounds on the integrals ( p.28|) - (|3.29| ) and ( p.35|) 
( [3.36| ). For example, for x large >> a, 



2 /■ /■ 1 

< x\GR\y > My)dy < -jMx) J <l)4y)My)dy = -'^M^) J g/_y_g/^y ^^- (3-45) 



From ( p.27|) we find, as x ^ cxd 



I s{-y - si+y'^y ^ ^'^ 



(3.46) 



I.e. 



oo 

/ < x\GR\y > (t>+{y)dy — > — 0+(a;). 
J gx 



(3.47) 



By using ( p.31|) and ( p.34| ), we see that as x -^ oo, apart from a relative constant normalization 
factor between ipeveniodd) and 0+, 



< x\GR\y > ipeven{y)dy — > (t>+{x)/gx 



(3.48) 



and also 



< x\GR\y > ipodd{y)dy — > (f)+{x)/gx. 



(3.49) 
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4. ITERATIVE PROCEDURE 

In this section, we present a systematic method to evaluate the eigenfunctions ipeven, i^odd 
and the eigenvalues Eeven and Eodd for the quartic potential case ( p..27| ). Write 

E = E -\- A 

-'-'even ^-' ^^ ^— *e 

and (4.1) 

Eodd = E + Aod, 

where E satisfies (|3.2| ) (and can be further tuned so that (|3.5|) holds). By definitions, ^even 
and ipodd satisfy 

\-'^ -^ JYeven ^eWeven 

and (4.2) 

{H - E)ipodd = ^od'ipodd- 



From ( [3.22|) , we have 

oo 

(H-E)- A, J < x\GR\y > ^,{y)dy = -A,^,{x), (4.3) 

— oo 

where j = even{e) or odd{od). Introduce 

oo 

Xj = ^j + ^j J < xIGrIv > ^j{y)dy, (4.4) 

— oo 

which, on account of ( [4-.2|) - ( |4-.3|) , satisfies 

{H-E)xj{x) = ^. (4.5) 

This is the same as equation ( ^31) for 0+(a;). Since both Xji^) and 0+(x) ^ 0, as x ^ oo, 
we have 

Xj{x)(X(p+{x). (4.6) 
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From ( p.48| ) - ( [3.49|) , the proportionality constant must not be zero, and can always be taken 



as unity by adjusting the normalization of ipj. Therefore, we can write 

oo 

1peven{x) = 4>+{x) - A^ J < x\G R\y > tpeven{y)dy (4.7) 

— oo 

and 

oo 

-ipoddix) = 0+(x) - Aod / < xIGrIv > ipodd{y)dy (4.8) 

— oo 

Impose the boundary conditions, at x = 

VLn(0)=0 and ^odd(0)=0 (4.9) 



From { ^.7\) - (|4.8| ) and regarding Gr as a square matrix, ^peven, ^odd, 0+ as column matrices, 
we derive 

i^eveniodd) = (1 + Ae(od)G r)~'^(J)+; (4-10) 



I.e. 



V^e.en = (1 " A,G R + A^GJ, - ■ ■ 00+ (4.11) 



^Podd = (1 - AodGR + Ai^Gj, - ■ ■ 00+ (4.12) 



Combining (^) with (|4ll| ) and using 0'+(O) = -0'_(O), we find 





,2A, 



oo 

2A r 

1 + ^ l[<Pl{y) + My)<P-iy)]dy 



*e\2 



[<t>+{yi) + <P-{yi)]dyi / [0-(i/i)0+(i/2) - (t)+{yi)(p-{y2)](p+{y2)dy2 



yi 

oo oo oo 

2 



-oo — oo 



+ T^l J [4>+{yi) + (t>-iyi)]dyi J dy2 J dys < yi\GR\y2 >< y2\GR\y3 > 0+(2/3) 

+ ••■ =0. (4.13) 
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Likewise, combining (|4.9| ) with ( |4.12|) , we derive 



1-^^ l[^4y)-^_{y)]^4y)dy 



A 







.00 00 

+i.^^f j[<P+i.yi) - <P-i.yi)]dyi j[(t^~i.yi)<P+i.y2) - <P+{yi)<P-{y2)]<P+{.y2)dy2 

S/l 



00 00 



-j^ld j [<P+{yi) - <P-{yi)]dyi J dy2 J dys < yi\GR\y2 >< y2\GR\y3 > (p+iys) 

+ ••• =0 (4.14) 



-00 —00 



From ( |4.13|) - (|4.14| ), Ae and Aod can be evaluated to successive orders of e". 
For a systematic analysis, we expand 

2A, 



A 



5o(e) + e5i(e) + e'62{e) + ■■■ + e"5,(e) + ■ ■ 



2A 



od 



A 



Define 



5o{od) + eSi{od) + t^52{od) H h e''5„(orf) + 



A 



< x\fR\y > = - < x\GR\y > 

0_(x)0+(y) -0+(x)0_(|/), x<y 
0. X > y. 

Cast (|4ll| ) and ( F34| ) into the forms 

1 + ao(^) + 6ai(^)^ + eV(^)^ + ■ ■ ■ + 6"a.(^)"+i + ■ ■ 



and 



l + 0o(— ^) + eOi 



A 



A 



■"^^2 + ^262 



,2A 



odx3 



A 



+ 



9 A 



where 



ao = y [0+(l/) + (t)^{y)(t)+{y)\dy, 




(4.15) 



(4.16) 



(4.17) 



+ ■■■ = (4.18) 



(4.19) 
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bo • '-' 



[<Pi{y)-<P^{y)<P+{y)]dy, 



(4.20) 



oo oo 

eai = j [(p+iyi) + (t)-{yi)]dyi j <yi\- fR\y2 > (p+{y2)dy2, 

-oo 



(4.21) 



e6i = - j[(t)+{yi) - 0_(i/i)]rfyi j <yi\- fR\y2 > (p+iy2)dy2, (4.22) 



e-a„ = / [(/)+(yi) + (j)-{yi)]dyi / < yi\{-fRT\y2 > 4>+{y2)dy2, 





(4.23) 



e'^bn = - J [Mvi) - (f>-iyi)]dyi I <yi\i-fRT\y2> My2)dy2- 





Theorem 



a„ = 0(l) and 6„ = 0(1). 



Proof. Write 



e"a„ = (-) V(+) and e"6„ = -(-)"J(-), 



where, in accordance with ( |4.23|) - ( [^.241) , 

oo oo 

J(±) = dyi dy2--- / dy^ / rfyn+i-^(±) 
and 



oo oo oo oo 



2/1 Vn-l Vn 



/(±) = [0+(yi) ± 0„(yi)] ■ [0_(yi)0+(y2) - Myi)4>-{y2)] 

■[(f)-{ym)(p+{ym+l) - 0+(ym)0_(ym+l)] " " " 
■[<P-{yn)<P+iyn+l) - 0+(yn)0-(yn+l)] " 0+(y„+l). 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



(4.28) 
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Decompose the integration range of ( |4.27|) into the sum 



(0) + (1) + (2) + ■ • ■ + (m) + ■ ■ ■ (n) + (ri + 1), 



where 



(0) : a < yi < y2 < ■ ■ ■ < Vn < Vn+i 



(4.29) 



(m) : <yi <y2< ■■■ <ym<Oi, 



a < ym+i < ■■■ <yn< 1/n+l 



with m = 1, 2, ■ ■ ■ , n, and 



5 -"5 5 



(ra + 1) : < 2/1 < 1/2 < ■ ■ ■ < 2/n < 2/n+i < « 



The integral ( |4.27|) can be written as 



n+l 

J(±) = Y. Jm{±). 
m=0 



where 



n+l 
'(m) ^ 



a a 



o oo 



(4.30) 



(4.31) 



(4.32) 



Jrn{±) = J n dyj{±) = j dyi j dy2--- j dy^ j dy^+i j dym+2 ■■■ J dy^+i ■ /(±). (4.33) 



J/l y-m-l OL Vm + l Vn 

Consider first Jq. From ( [3.44| ) and (|4.28|) , the magnitude of /(±) satisfies the inequality 



|/(±)|,„(o) < i-r\Myi) ± 0-(i/i)l[0+(i/i)04y2)][0+(2/2)0-(i/3)] ■ • • 

[0+(Z/n)0-(yn+l)]0+(2/n+l). (4.34) 



Since at any x, on account of ( [3.25| ), 



(xUJx) = 0i 



(4.35) 



and for x > —a, 
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0^(x) < 0(1), (4.36) 
we have 

|0+(l/i)±0_(2/i)l0+(2/i)<O(l) (4.37) 
and therefore 

J^=o(±) = 0(e"). (4.38) 

Next, for any m > 0, in the region (m) given by ( [^.301 ), we have, on account of (|3.41| ) and 

\I{±)U^) < (-)'^-'"[0+(yi) ± <P-{yi)] ■ \<P-{yi)My2)] ■ [0-(y2)0+(y3)] ■ ■ • 

[(l)-{ym-l)<P+{ym)\ ■ \(t)-{ym)<P+{ym+l) " 0+(ym)0_ (z/m+l) | " 
[0+(ym+l)0-(ym+2)] ■ ■ • [0+(yn)0-(yn+l)]0+(2/n+l)- (4.39) 

For < yi < a, we have 0+(i/i) > 4>-{yi) > 0, 

Myi)<PAyi) = oie) 

and therefore 

[Myi)±^-iyiM-iyi) = oie). 

Along the sequence on the right hand side of ( [4.39|) , we find 

My2)<PAy2) = oie),--- 

^+{ym-i)<P-{ym-i) = 0(e), 
and since both (j)+{yni)4>-{ym) and (f)+{ym-i)4>-{yrn-i) are 0(e), 

28 



0+(Z/m)[0-(l/m)0+(Z/m+l) " 0+(ym)0_ (l/m+l)]0+(l/m+l) 
= 0+(Z/m)0-(l/m)0+(l/m+l) - 0+(ym)0- (y,n+l)0+(ym+l) 

= 0(e). 
Likewise, 

0-(ym+2)0+(2/m+2) = 0(e), 

(f)-{yn+i)(f)+{yn+i) = 0{e). 

Thus, in the region (m), /(±) = 0(e"), so is Jm(±); i.e., 

J^(±)=0(e"). (4.40) 

Likewise, we can show J„+i(±) = 0(e") and that leads to (|4.25| ). QED 

Substituting ( |4.15| ) into ( |4.17|) , and setting the coefficients of different powers of e sepa- 
rately to be zero, we have 

l + ao5o{e)=0, (4.41) 

ao6i{e) + aMef = 0, (4.42) 

aoS2{e) + 2aiSo{e)5i{e) + a2(5o(e) = 0, etc.. (4.43) 

Likewise, 

1 + Mo(orf) = 0, (4.44) 

Mi(orf) + bi6o{od)^ = 0, (4.45) 

boS2{od) + 2bi6o{od)5i{od) + 02^0(0^) = 0, etc.. (4.46) 
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(We recognize that, in view of ( |4.19|) - (|4.24| ), a„ and 6„ also contain e-dependence,, although 
each is 0(1) to the leading order. We likewise permit the (5„ to have implicit e-dependences so 
that they may not be fixed uniquely by (|4.15| ). We have removed the ambiguity by requiring 
that when ( |4.15|) is substituted into ( |4.17|) - ( |4.18|) , the resulting equation should be true to 
each explicit power of e, disregarding the implicit entrance of e through the a„, 6„ and Sn- 
The important point is that if Ag, A^^ are given by ( [4.15| ) with ( |4.41| ) - (|4.46| ), then (|4.17| ) - 
( p^ ) will follow.) 
To 0(e), 



Since 



and 



2Ae = 


= A5o(e) 


= - 


-X/ao 


2A,, = 


= X5o{od) 


= - 


-A/60. 


oo 



oo 


Oie) 





We find, from ( |4l^ ) - (^[^g). 



M+{y)dy 



ao= 4>l{y)dy + 0{e) 



(4.47) 



(4.48) 



4>l{y)dy + 0{e) 



Thus 



A, 



00 
rA / f - 2 



and 



which lead to, neglecting 0(e^ 



^J !<Pl{y)dy] + 0{ 



Ao. = +[^J j<t>l{y)dy] + 0{e^ 




-{Eeven + Eodd) — 90. — "T^ 



9 



89 



4a2 2g{2af 2g^{2a) 



+ o{-j-^) 



g^a' 



-A. 



A, 



A 



od 



X = Sga^e—s^" [1 



53 



8ga^ 256(^a3)2 
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+ 0{ 



{ga^y 



dx. 



01 



(4.49) 



(4.50) 



and 0+ = e-^(+) given by (^) and ( pl ). 

By using the asymptotic expansion of S{+) in 0+ = 6""^*^+^ 

/dv 

and ( [^.41] ) - (|4.46| ), we can evaluate the coefficients of (ga^y^'e^'^^^'^ in the double series 
expansion ( |1.33|) - (|1.34| ) for Eeven and Eodd in terms of definite integrals. In particular, the 
corresponding asymptotic expansion ( |2.12|) of E is automatzca//?/ agreeing with ^(Eeven+Eodd) 
to 0(e). The underlying reason is: from ( |3.13|) , a variation of a —^ a + Sa gives Scp^/Sa < 
0(e), which leads to < O(e^) to Ag and Aod 

As remarked in Section 3, different choices of ±a within the range (|T^ ) give different E 
and 0-1- within O(e^); however, they lead to the same i>even(odd) and Eeven(odd)- 
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5. REMARKS AND OTHER EXAMPLES 

(1) Ansatz of trial wave functions. 

In the notation of (pTl] ) - (p~4l), let all the minima of the positive potential 



V{x)=g^v{x) (5.1) 



be at points 



h,h,- ' ' Jn] (5.2) 



I.e. 



Vih) = 0, (5.3) 

with j = 1, ■ ■ ■ , A^. For each j, there exists a solution of the Hamilton-Jacobi equation ( p..9|) : 

i(^)2-.(x)=0 + . (5.4) 



According to ( p..l3|) , this solution can be written as 



X 

So{j) = ±J[2v{y) + {0+)]-2dy, (5.5) 

h 

with the ± sign chosen to make 

So{j) > everywhere. (5.6) 

These A^ 

^(j) = e-3^«(^') (5.7) 

functions are all linearly independent and can be used as a set of trial functions. In the 
neighborhood of x = Ij, the corresponding 0(j) describes the correct harmonic oscillator 
ground state behavior. The value of 0(j) at a different minimum x = k {i ^ j) gives the 
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barrier penetration amplitude from j to i, similar to the instanton description l^'^l. Therefore, 
the set (^]7|) forms a convenient ansatz of trial wave functions for the N low-lying eigenstates 
of the Hamiltonian (|1.2|). In the case when V{x) is periodic, from this ansatz we can construct 
the lowest energy band structure of the problem. 

(2) The systematic procedure developed in Section 4 can be extended to a much larger 
class of potentials with degenerate minima. The limitation depends on the applicability of 
equations similar to ( |4.7|) - ( |4.8|) . Here we give two examples that illustrate this limitation. 

Consider first the special potential 

V{x) = u[-6{x -I)- 6{x + /) + q6{x)], (5.8) 

where /, u and q are positive constants with 

q < 1, (5.9) 

so that 

oo 

f V{x)dx < -u. (5.10) 

The ground state of H = -\{<P/dx^) + V{x) is 

p l^\X I) T ^ J 

Yeveny^j ^ Yeven\ ^) ^ \ l^.iij 

(1 - «)e-«(^-0 + u^nix-i) 0<x <l 

^ K,' K ' 



with its eigenvalue 



It can be readily shown that 



Ee. = ~\^. (5.12) 



^ = ^[l + !^^!^e-2'^']. (5.13) 



K + qu 
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Thus, for q < 1, K \s > u, consistent with ( p. 101 ). If we neglect the barrier penetration factor 
g-2K/ ^i^^gi^^ K = u, as would be the case for a single —u5{x — I) potential. 
Introduce (p+^x) for this problem as the solution of 



:h-e)m^) = o, 



(5.14) 



with 0+(oo) = 0, 0+(/) = 1 and 



E 



— u 
2 



(5.15) 



We have 



-u{x—l) 



^u{x—l) 



0+(x) = < 



-u(x+l) -)- (1 — g)e"(^~'\ 

(g + 2(1 - g)e-2««)e"(^+^) + (-1 + g)e"2«'-«(^+0^ 
Construct the Green's function Gr in terms of 0+(x) and 0-(x) = 0+(— x) 



X > I 
< X <l 
-I < X <0 
X < -I. 



(5.16) 



< x\GR\y >= - 



where, similar to ( |^ ) 



-<P+{x)<P-{y) + <P-{x)<l)+{y) > for X <y 
, for X > y 



(5.17) 



A = 0'j_(x)0_(x) — 0'_(x)0+(x) = 2^(1 — q)e 



-2ul 



(5.1^ 



Consequently, as in ( [3.22| ) 



{H -E) < x\GR\y >= -5{x - y). 



(5.19) 



Because ii^e^ < -E, the integral 



< x\Gn\y > ipeven{y)dy 



(5.20) 
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is convergent. Therefore, as in ( |4.5|) , 

oo 
(H - E)[i)even + Ae J < x\G R\y > 1peven{y)dy] = 0, (5.21) 

— oo 

where Ae = E^v - E. However, unlike ( |47| ) - ( ^~8|) , we find that 

oo 

Xeven = i^even + Ag / < x\G R\y > 1peven{y)dy (5.22) 

— oo 

satisfies 

Xeven = 0. (5.23) 

For the quartic potential, instead of ( |5.23| ), we have (|4~6|), which can be written as 

oo 

Ipevenix) = 4>+{x) - A^ J < x\GR\y > 1peven{y)dy (5.24) 

— oo 

in accordance with (|4.7| ); furthermore 

i^even{x) ~^ (p+{x) 35 X ^ CX). (5.25) 

Recall that equations (|3.3C1| ) - (|3.32|) are applicable for any V{x). Assume 



9^ 
V{x) -^ —x"^ , as X ^ oo; (5.26) 



then instead of (|3.34| ), we have for large x 



For ra > 2, ^ is insensitive to E. The marginal case is ra = 2; in that case ( |5.24| ) remains 
valid, but both terms on its right hand side are important as x ^ oo. For the potential ( |5.8| ) 
of the above example, we have n = 0; therefore (|5.24| ) is not valid, being replaced by ( |5.22| ) 

(3) We will now examine the marginal case n = 2. Assume 
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V{x) = g^v{x) = < 



Ig'^ix - ly for X > 

k5{x) 
\g^{x + iy for X < 0. 



Let 6, be the solution of 



ff0+ = --0'|(x) + V(x)0+(x) = E0+(x) 



Write 



(5.2^ 



(5.29) 



[X] = e 



-si+) 



(5.30) 



with 0+(oo) = 0, and expand, as in ( p..6|) - ( p..7|) , 



5(+)=^5o(+) + 5i(+)+^-^52 



(5.31) 



E = (7^0 + El + (7-^E2 + ■ ■ 



(5.32) 



From ( [5.29|) , we have 



S{+)" + S{+)" + 2g'v{x) = 2E. 



(5.33) 



Substituting ( |5.31|) and (|5.32| ) into (|5.33| ), and setting the coefficients of different powers of 
g to be zero, we derive 

Soi+y = ix-i), 

which, together with the "normalization" condition 0+(/) = 1, gives 
in addition. 



-^0 = ^ ; El = E2 = ■ ■ ■ = 
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and 



S^i+) = S2i+) = ■ ■ ■ = 0. 



Thus, 



(x) = e-59(^-')' for X > 0. (5.34) 



At X = 0+, 



.l„l2 



and (5.35) 

0'^(O+) = gle--2^^\ 



On account of the A5{x) term in V{x), at a; = 0— we have 



^9l' 



</.+ (0-) = e-2^ 

and (5.36) 

(/.V(O-) = -2A0+(O)+0',(O+) = (-2A + (7/)e-^5^^ 



For X < 0, (f)+{x) satisfies the second order differential equation ( |5.29D 

-^<f>l + ^g'{x + in+ = ^g<P^, (5.37) 

whose general solution is 

0+(x) = ax{x) + 6e~59(^+0^ (5.38) 

where a and b are constants and 

a{x+/)2 l'g{y+lf 



X{x) = e-t(-+')^ / e3^y+'> dy. (5.39) 



Because of the boundary conditions (|5.36| ), 

a = 2[gl - A]e~^'' and 6=1. (5.40) 

Hence, for x < 
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[x] = e 2 



hai^+if 



[1 + a fe'^y+^^'dy]. 



(5.41) 



Introduce, as before, 



0_(x) =0+(-x). 



(5.42) 



we have, for x > 0, 

X 

(j)_{x) = e-5s(-0^[i _ a fe^'-y-^^'dy]. 



When A = gl, we have a = and the exact groundstate wave function given by 

■ig{a;-/)2 



(5.43) 



e 2 

_i 
e 2 



X > 



^9(x+lf ^ ^ Q 



the corresponding groundstate energy is 



2'- 



For A 7^ (jf/, the Hamiltonian H can be written as 



(5.44) 



(5.45) 



H = Ho + h, 



(5.46) 



where 



and 



Hn 



"2dx2 



+ < 



gl5{x) 



bH^ + i?, 



X > 



X <0, 



(5.47) 



h = {A-gl)6{x). 



(5.48) 



Regarding h as the perturbation, to the first order in h, the change of energy is given by the 
expectation value of h, with (|5.44| ) as the unperturbed wave function; i.e.. 
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Ee.en = l9+^4^^^^^. (5.49) 

^ Of p-g{x-lYdx 



We now turn to an alternative derivation of E even by using (|4.7| ) and ( |4.9|) . Let (p+ix) and 
0_(a;) be the same functions that satisfy ( |5.29| )-( p^ ). From (|5.34| ) and (|5.43| ), it follows 
that for X and y both positive definite, 

0+(x)0_(y) - <P-{x)My) = aMx)My)[fix) - fiv)], (5.50) 

where 

X 

f{x)= fe'^'-^^'dz, (5.51) 



and 

<f>'4x)<j>4y) - <jyL{x)My) = -9{x - i)[M^)<f>4y) - <l^4^)My)] + ^^44 (^-52) 

0+(x) 

which, at X = y, becomes 

0'^(x)0_(x) — </)'_(x)0+(x) = a. (5.53) 

As in (|5.24| ), we write 

CO 

Ipevenix) = (j)+{x) - Ag / < x\GR\y > 1peven{y)dy, (5.54) 

— oo 

where 

Ee.en = f + ^e (5.55) 



and 



, l{^~ix)4>+{y) - ^+{x)^-{y)] , x<y 

< x\GR\y >= { (5.56) 

, X > y. 
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For y > X, hy using ( |5.52|) - ( |5.53| ) we see that as x — > 0+ 



£ < x\Gn\y >^ gl < 0|G«|y > -2[0+(|/)/0+(O+)]. (5.57) 

The function ipevenix) and 0+(x) are continuous at x = 0, but not their derivatives. We have 

0^(0) = e-^^'' (5.58) 

oo 
IpeveniO) = C'^^^^^ - A^ J < 0|Gij|y > ^even{y)dy (5.59) 

— oo 

0'^(O+) =^/e-55'' (5.60) 

and 



oo 



Cen(0 + ) = gl^eveniS^) + 2Aee29' j 0+(x)^e.en(x)cia;. (5.61) 



From the Schroedinger equation with V{x) given by ( |5.28D , the groundstate wave function 
must satisfy 

^Len(0+) V^Len(O-) 



Substituting ( |5.61J ) into ( |5.62| ), we derive 



A. (5.62) 



oo 



O A 

^-91 = TTTre^^'^ / (j)+{x)ilJeven{x)dx. (5.63) 

To derive Ag to 0(e~^' ), we need only the zero*^ order approximation in the wave function 

^even = 0+ (5.64) 

which gives 



^^^ (A-,/)e-»" ^ , 



'e - 00 I -V- ; (5-65) 

2/02_(a;)da; 




in agreement with (|5.49| ). 
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To derive the next order correction, we replace ( |5.64| ) by 

oo 

'ipevenix) = (p+ix) - Ag / < x\GR\y > (j)+{y)dy. (5.66) 

— oo 

From (|30D, (|33D and (|536| ), we have for y > x > 0, 

y 
< x\GR\y >= 20+(x)0+(y) J ^-\z)dz. (5.67) 

X 

Hence, the approximation (|5.66| ) leads to 

oo y 

1peven{Q) = (t>+mi - 2/\e j 4>l{y)dy j r+\z)dzi (5.68) 



oo y 

i^:LiO) = <P^\0)[l + 2A,l<Pliy)dyJ<p-'iz)dz] (5.69) 



and 

oo oo oo 

jM-)i:U-)d- ^Jd.M.)lM-) - A./ < xiG.l, > My)dy] 

X 

oo 

oo j (j)\{x)F{x)dx 

= f dy<l>l{y)[l - 2A,^ ], (5.70) 

I4>Uz)dz 



where 

oo y 

F{x)= j<Pl{y)dyjr^\z)dz. (5.71) 

X X 

Thus, ( |5.69| ) can also be written as 

C.en(O) = 0;'(O)[l + 2AeF(O)] (5.72) 

which, together with ( |5.70|) , gives 

oo oo 

(t)+{x)%l)even{x)dx = / 4>X{y)dy 



-ipeveniO) ^ 0+(O) ^ 

2A, 



7 \ r 

■{1 - ^ ^— / <Pl{x)[F{x) - Fmdx}. (5.73) 
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Because the potential V{x) -^ \g^x'^ as a; ^ oo, which is the marginal case ra = 2 in 
accordance with ( |5.27| ), the corresponding integral F{x) given by ( |5.71J ) has a logarithmic 
divergence as the range of the y-integration -^ oo. However, as we shall see, the difference 
F{x) — -F(O) and, therefore, Ag are well-defined. 

We modify ( [5.71|) by introducing a regulator (1 + \y)^^ in its integrand. Let 



Fa(x) = / (1 + \y)-^<Pi{y)dy / <p-^\z)dz (5.74) 

and define 



F{x) - F(0) = lim \Fx{x) - Fa(0)]. (5.75) 



In ( |5.74D , the 2;-integration is 



y y 

2 



{z)dz = / e^^""') dz 



X X 

2 



(^9{y-i? Qg(^-if } q9(z-i) 



+ / 7— ^c?z. (5.76) 



2giy-l) 2g{x-l) 7 2g{z - If 

The subsequent y-integration for the first term on the right-hand side of (|5.76| ) leads to: 

p9{y-i? 
'l + \y)-^ct?,{y)— -dy 



'igiy - /) 

^ -[InX + In^^^^]. (5.77) 



2^(1 -I- AO l + Xx' 

In the difference Fx{x) — -Fa(O), the above InX term cancels. The limit A -^ 0+ gives a finite 
F(x) — F{0), therefore, a finite ratio 

1 



ll^even{0) 



(l)+{x)ljJeven{x)dx (5.78) 





in accordance with (|5.73| ). The substitution of ( |5.73|) into (|5.63| ) leads to a well-defined 



quadratic equation for Agi its solution determines Ag to the 0(e ^^' 
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APPENDIX 

In this Appendix, we discuss further the difference between a ^/'-function and a 0-function, 
introduced in Section 1 (between (|1.8| ) and (|1.9|). For clarity, we return to the quartic potential 
case in which Hamiltonian H is given by (|1.27| ) and its even and odd groundstate wave 
functions are iJeven and 'ipodd- As in ( |1.29|) , 

Mx) = -[tpevenix) ± i)odd{x)]. (Al) 

Let (f)±{x) and E satisfy ( p.l| ) - ( p.6|) . At x = ±00, ^±(x) are well-behaved, with 

V^±(oo) = V±(-oo) = 0. {A2) 

However, since 

{H - E)(f)±{x) = 0, {A3) 

with E ^ an eigenvalue of H, only 

0+(oo) = 0_(-oo) = O, (AA) 

but 0+(— 00) and 0-(oo) are divergent. 

From ( |4.1j ) - (P~5|), we express (P~^), Xj(x) oc 0+(a;), as 

Xj(x) = 90+ (x), (A5) 



where Cj are constants and, as before, j = even{e) or odd{od). Thus, ( p~7|) and ( |4~8|) now 
become 



'ipevenix) = Ce0+(x) - Ae / < xlG/j]?/ > ll)even{y)dy {A<o) 

—00 

and 
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-ipoddix) = Cod(l)+{x) - Aod / < xIGrIv > ipodd{y)dy. {A7) 

— oo 

Through x -^ —x, it follows that 

oo 
i^evenix) = Ce(j)-{x) - Ag / < x\GL\y > 1peven{y)dy (AS) 

— oo 

and 

oo 

i'oddix) = -Cod<l)-{x) - ^od J < x\GL\y > ipodd{y)dy. {A9) 

— oo 

Combining (Al) with (A6) - (A7) and using the matrix notation of ( |4.1C]| ), we derive 



and 



ip. = c-cf>+ - G'r[A_^+ + A+^_] (All) 



where 



C± = |(Ce ±Cod) 

and (A12) 

A± = i(Ae±A,rf). 

Likewise, from (A8) - (A9) 

V'+ = c_0_-G'l[A+V'+ + A_V-] (^13) 

and 

^_ = c+0_ - Gl[A_^+ + A+^„]. (A14) 

On account of (|3.48| ) - ( p.49|) and the corresponding formulas, replacing Gr by Gl, we find 
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as x ^ oo, 'ip+ix) -^ c+0+(a;) -^ (^15) 

^_(a;) -^ c_0+(a;) ^ (A16) 



and 



as a; -^ — oo, ■?/'+(a;) ^ c_0_(a;) ^ (^17) 

^_(a;) -^ c+0_(a;) ^ 0, (A18) 

in agreement with (A2). As we shall see, it is convenient to choose 



c+ = 1 and c_ = e = e's^" . (A19) 



Thus, 



d+ix) as a; — > oo 

^+(x) ^ <( ^^^ ^ (A20) 

e0_(x) as X ^ — oo 



and 



ech+ix) as a; — !► oo 

^_(a:)-.<( "^^^ (A21) 

(j)-{x) as a; — !► — cxd. 

We now examine the expressions 

0^ = e-^(±) (A22) 

and the asymptotic expansions (neglecting 0(e) corrections) 

S{x) = gSo{x) + Si(a;) + g''^S2ix) H 

^(x) = gSoix) + 5i(x) + g-'S2{x) + ■■■. (A23) 

From (A20) and (A21), we see that 
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S{+) 



4„^3 



ga' + Si-) 



as a: — i> oo 
as a; ^ — oo. 



(A24) 



and 



as a; ^ oo 
as a; — s> — oo. 



(A25) 



Iga' + Si+) 
s{-) 

Since, neglecting 0(e), {S„(±)} and {^^(i)} satisfy the same set of differential equations 
we have 



^ra( + j — Jn[ 



and 



4 
— ( 
3 



Soi+) = ^a^ + SoH 



^m\') ^m\ 



Likewise, 



^n{ } — ^n[ 



and 



4 
— ( 
3 



Soi-) = ^a^ + Soi 



^m( ) — Jm[~r] 



for x>-a + 0{^), (A26) 



ga' 



for X < -a {A27) 

for X < —a and m 7^ 0. {A2^ 



for a;<a-0(-i^), (A29) 



for x> a (A30) 

for X > a and m 7^ 0. (^31) 



We note that So(±) and S'o(±) are both solutions of 



2 „2\2 



[S[,(±)]^ = 2t; = (a;'-a' 



(A32) 



and 



2 „2\2 



[5;(±)]^ = 2t; = (a:^-a^) 



{A33) 
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In accordance with ( |2.11|) , at all 



X 



Sq{+) = l(x-a)2(x + 2a) 



and 



(A34) 



So(-) = Ux + a)H-x + 2a) 



However, from (A26) - (A27) and (A29) - (A30), 



So(+) 



and 



So{+) = \{x-ay{x + 2a) 



for X > —a 



|a^ + S'o(-) = |a^ + ^(x + a)^(-a; + 2a) for x < -a. 



fa^ + S'o(+) = |a^ + \{x - a)2(x + 2a) 
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|(x + a)^(— X + 2a) 



for X > a 
for X < a. 



which are solutions of the Hamilton-Jacobi equation 



(A37) 



(v435) 



(A36) 



with 0+ as the total energy. By requiring So(±) = oo at x = ±cx) and So(±) = at x = ±a 
respectively, we derive (A35) - (A36) in the limit OH > 0. 
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